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Abstract 

Using confocal microscopy, we measure the packing geometry and inter- 
droplet forces as a function of the osmotic pressure in a 3D emulsion system. 
We find that the nonhnear elastic response of the pressure with density is not 
a result of the anharmonicity in the interaction potential, but of the correc- 
tions to the scaling laws of the microstructure away from the critical point. 
The bulk modulus depends on the excess contacts created under compres- 
sion, which leads to the correction exponent a = 1.5. Microscopically, the 
nonlinearities manifest themselves as a narrowing of the distribution of the 
pressure per particle as a function of the global pressure. 

The elasticity of jammed sphere packings sheds light on the broader question of 
how amorphous solids support stress ifTTlfTSl . In the limit of purely repulsive fric- 
tionless particles, applicable to emulsions and foams, theoretical arguments sug- 
gest a singularity at the critical random close packing density <^c [27, 29 1. This 
singularity gives rise to the surprising scaling behavior of the elastic moduli and 
the microstructure, as observed in numerical simulations lH [171 HOI 121 IS l- In 
particular, the excess particle contacts created upon compression are thought to 
control the singularity in elasticity. However, an experimental test of the scaling 
law between the coordination number and the applied pressure is lacking in the lit- 
erature. This is because it is experimentally difficult to access the regime very close 
to the jamming transition with enough resolution to test the scaling laws. More- 
over, the contact network in 3D is typically hidden from view. On the other hand, 
bulk measurements have shown that compressed emulsions exhibit deviations from 
the predicted linear scaling of the osmotic pressure with the global density above 
random close packing [18]. This scaling was conjectured to arise from the contact 
number dependent anharmonicity observed in the potential of a single compressed 
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droplet lfT4l . While a stiffer potential could lead to nonlinear elasticity in the limit 
of marginal rigidity ||29]| . the making and breaking of contacts are other plausible 
sources for the deviation f24l. A microscopic experimental approach is therefore 
needed to decipher the origin of the macroscopic elastic response of the system. 

Here we confocally image compressed emulsions to measure the packing ge- 
ometry and the interdroplet forces in 3D, from which we then infer the scaling 
laws between the applied pressure P, the average number of contacts z, and the 
global density (|) away from the critical point. The statistical fluctuations of these 
quantities within the packing give us additional insight into the mechanisms of 
stress transmission in this amorphous medium. The oil-in-water emulsion sys- 
tem described in references lITSl 13 is athermal with an average droplet radius 
{R) = 2.5 ^im and a polydispersity of 25% preventing crystallization [22|. The 
refractive index matched emulsion is transparent and the Nile Red dye reveals the 
3D packing of droplets using a confocal microscope (Leica TCS SP5 II). The low- 
est compression rate is achieved by creaming under gravity, which gives a packing 
fraction of = 0.68 ±0.02, while centrifugation for 20 minutes at an acceleration 
rate of 3000g leads to foam-like structures with (|) = 0.88 ±0.02, as shown in Figs. 




Figure 1: Confocal 3D images for emulsions creamed under gravity in (a) and 
compressed at P = 2.2kPa by uniaxial centrifugation in (b). (c) A zoom of two 
reconstructed droplets labeled a and P shows their geometric overlap area A"P, 
which is used to estimate the interdroplet force. 
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Figure 2: Scaling laws for the excess average number of contacts 5z with the 
change in global density 8(|) = ((|) — (|)c ) in (a) and the applied pressure 5P= (P — Pc) 
in (b). A parametric plot of 5P versus 5(|) in (c) implies the nonlinear scaling law, 
in agreement with previous bulk measurements. The power-law fits are shown on 
log-log scales in the insets, (d) Probability density distributions of the pressure per 
particle rescaled by the mean as a function of the applied pressure. The inset 
shows the corresponding coefficients of variation, CV for all measured distribu- 
tions. The pressure fluctuations aie best fit with exponential tails at low pressures 
and a Gaussian distribution at high pressure. 

[l]\,B, respectively. Allowing the emulsion to relax to its uncompressed state over 
a period of several days probes a broad range of intermediate packing densities. 
The packings are analyzed using a Fourier transform algorithm to reconstruct the 
packing ||4l and a navigation map tessellation to identify which particle resides in 
each cell volume l[2ll6l[T2l. 
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For each particle, we next obtain the local volume fraction (j)L, the number 
of contacts ^^'^ the pressure on that particle P^. The statistics of the local 
parameters in a given packing yield the corresponding global parameters. The 
local volume fraction is defined as the ratio of the volume of the particle to the 
volume of its navigation map cell, while the global packing fraction is the total 
volume of the droplets divided by the volume of the box. To determine the contact 
network, we test whether there is a geometric overlap between the reconstructed 
spheres (5), as shown in Fig. \Vp. This procedure leads to an estimate of the mean 
contact number z it 0.3, where the uncertainty is due to the resolution limit of the 
microscope. Moreover, the area of overlap, A"P, is proportional to the repulsive 
interdroplet force F"^, which is given by the Princen model 11231 : 

77 «P = (1) 

where the weighted radius /?"P = and the interfacial tension a = 9.2 mN/m. 

For small deformations, this model gives an essentially linear force law. This re- 
sult is corroborated by the fact that the geometric overlap A"P agrees with the area 
of highlighted fluorescence intensity between the droplets, i.e. the area of defor- 
mation [3|. Moreover, minimizing the energy of this system of linear springs to 
satisfy force balance does not move any of the particles beyond the resolution of 
the technique, which is one third of the voxel size, i.e. lOOnm. Note that the in- 
compressibility of the emulsion leads to a divergence of the pressure at very high 
compression. In addition, the ability to lose or create contacts is another source of 
nonlinearity ll24l . 

The measurement of the forces then allows us to calculate the local pressure per 
particle P^, which is equal to the trace of the corresponding Cauchy stress tensor, 

where indices i,j represent vector projections of each contact force F"P, is the 
volume of the navigation cell of the droplet a, r"P is a vector from the center of a 
to its point of contact with the droplet P, and the sum of the force moment tensor is 
taken over all contacting droplets. Applying Eq. (2) to the entire packing gives the 
global pressure P. We find that the diagonal elements of the global stress tensor are 
close to each other in magnitude, suggesting that the compression of the emulsion 
is isotropic due to the absence of friction. We repeat these measurements on 50 
samples with global pressure values ranging from 0.27 ±0.05 kPa under gravity to 
2.29 ± 0.05 kPa. 

Let us first assess how the excess contacts grow with the global density. Our 
system of frictionless spheres reaches isostatic equilibrium at the critical point. 
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Figure 3: Probability density distributions of local parameters (^) ^^id '^h 
(b) are rescaled by the mean for three different global pressure values. Insets show 
the corresponding coefficients of variation, CV, for all measured distributions. The 
data are weU fit with the granocentric model. 

which fixes its average coordination number to zc = 2<i, where J = 3 is the dimen- 
sion [TII29I. To deduce the volume fraction we take into account the shape defor- 
mation of each droplet away from its spherical shape, such that the droplet volume 
cannot be larger than that of its cell. In Fig. 2A we show the scaling of excess 
contacts with the square root of the packing density 5z = Z — Zc = Zo y/W^^c) = 
Zo \/5^ over the experimental range. This geometric result is in good agreement 
with previous numerical simulations and theoretical arguments |ill|2T]|291, as well 
as 2D data on frictional particles |16 | and foams |13|. Moreover, the value of the 
prefactor zo = 10.6 is in good agreement with the one obtained from numerical 
simulations of bidisperse spheres close to the jamming point, zo = 8.4 it 0.5 [21]. 
This fit determines a higher value of (|)c = 0.68 for the polydisperse packing than 
the predicted monodisperse random close packing density, (|)c = 0.64 [25|. The 
obtained value is in excellent agreement with the density measured at the lowest 
pressure. 

By contrast, the scaling of excess contacts with the applied pressure in Fig. 2B 
cannot be fit by the square root law predicted for a network of harmonic spheres 
close to the jamming point |29J: 5z ~ VP — Pc, where Pc = 0.27kPa is the residual 
pressure under gravity. This scaling behavior was observed earlier in simulations 
of harmonic, frictionless particles HllllJ. Surprisingly, such a scaling was ob- 
served experimentally in the case of frictional disks over a narrow range of 1% in 
density 1.16.1 . Our experiments are the first to test this relation in frictionless pack- 
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ings and in 3D, over a much broader range of packing densities up to 0.2 above (|)c. 
The data shows important deviations from the scaling prediction, explained below. 

Combining the two square root laws for the scaling of excess contacts with 
density and pressure predicts a linear parametric relationship between pressure and 
density, where 5P = P — Pc = Po5(|)^ \ is the prefactor and y = 2 is the expo- 
nent in the harmonic interaction potential |[2Tll27l . However, the deviations of the 
data from the square root law in Fig. 2B lead to the nonlinear relationship of pres- 
sure versus density presented in Fig. 2C. These measurements are in good agree- 
ment with those obtained from bulk experiments ifTSl . which measure the pressure 
macroscopically. The agreement between the two data sets gives additional validity 
to the harmonic assumption in Eq. (1) used to calculate the pressure in our system. 
Therefore, the observed deviations from the linear scaling are likely to come form 
the formation of contacts upon compression and not the proposed anharmonicity 
in the potential 114]. 

Theoretically, these rearrangements can be understood in terms of the correc- 
tions to scaling away from the critical point. It has been shown that for generic 
random elastic networks, all the elastic moduli must vanish linearly with the ex- 
cess coordination G ~ B ~ 8z |[27ll28l . as is observed numerically 191. However, 
packings of particles differ from random networks in the following sense: their 
geometry is constrained by the fact that all the contact forces between the par- 
ticles are positive. It was shown that this constraint implies that the bulk mod- 
ulus must have an additional contribution that does not vanish at jamming, that 
is B = (|)3P/3(|) Ci +C25z, where Ci andCi are constants. Moreover, the pa- 
rameter 6z has been shown to govern the crossover from the isostatic behavior 
close to jamming to the continuum behavior at the large scale |[T0|| . Since the 
scaling relation 5z ~ \/5^ holds over the experimental range in Fig. 2A, we get 
B ^ Ci +C'2\/^ [28 1 . If we assume that, apart from the creation of new con- 
tacts, no plastic rearrangements occur, this equation applies for all (|) and can be 
integrated and expanded to give the two leading order terms in 6(|), 

P-Pc=Poi<\>-<\>c)+Pi{<\>-<\>cT (3) 

where Pc = 0.27 kPa is the experimental value of the pressure at the lowest com- 
pression, Pq = 3.0 kPa and Pi = 14.8 kPa are prefactors obtained from the fit of the 
data in Fig. 2C, and a = 1.5 is fixed as the universal exponent for the first correc- 
tion to scaling. The data is in excellent agreement with the theory in the range up to 
S(|) = 0.18, beyond which the incompressibility of the droplets causes a divergence 
in the pressure. Remarkably, converting the prefactor Pq in the linear term to the 
units of pressure defined in numerical simulations of bidisperse packings close to 
jamming recovers the same value to one decimal place 1,21,1 . The success of the fit 
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in Fig. 2C with Eq. (3) implies that the data in Fig. 2B can also be explained in 
terms of the corrections to the scaling of excess contacts with pressure in the form 
8z = Di\/hP-D2hP, where the prefactors Di = 4.3 kPa"'/^ ^nd D2 = 0.8 kPa"^ 
are obtained from the fit in the figure. The prefactor values agree with the para- 
metric prediction of the fits in Figs. 2A,C, which testifies to the accuracy of the 
experimental data and the validity of the model. 

Microscopically, the distributions of the pressure per particle for three dif- 
ferent pressures are shown in Figure 2D. The tails of the distributions narrow down 
from exponential to Gaussian as a function of the applied pressure, as shown in 
Fig. 3C. This is consistent with theoretical considerations for the distribution of 
forces close to the jamming point Q 111 |26l and the randomization of stress that 
occurs as the pressure is increased |l30l. The change in shape of the local distri- 
butions with the applied pressure contributes to the nonlinear elastic response of 
the material. Finally, we show the microscopic distributions of (Sf^ as a 

function of pressure in Fig. 3. It is interesting to note that rescaling the contact 
distribution by the mean collapses the distributions at different pressures. This 
is also confirmed by the invariant coefficient of variation shown in the inset. By 
contrast, the distributions of local packing fraction do not collapse with the pres- 
sure, but become more peaked. The distributions in Figs. 3A,B are superimposed 
with the theoretical predictions of the granocentric model [6','T9l, which takes the 
measured average number of contacts, neighbors and the global density as inputs 
and generates the local fluctuations in each input parameter. This model is based 
on two local random processes: filling the space around each particle with solid 
angles contributed by the neighbors and choosing some of those neighbors to be 
contacts to ensure rigidity. The good agreement between the data and the model 
suggests that random processes on the level of the first neighbor shell are sufficient 
in explaining the diversity of local configurations. 

In conclusion, our data show corrections to the linear scaling of the pressure 
with density that arise because the bulk modulus is not constant near jamming, but 
depends linearly on the excess contacts in the force network. We also show that the 
interaction potential is nearly harmonic, independent of the compression. There- 
fore, the first correction to scaling, which only takes into account the creation of 
contacts but not particle reanangements or the incompressibility of the droplets, is 
sufficient in accounting for the nonlinear dependence over a wide range of densi- 
ties. This study reveals the mechanism by which emulsions are stiff to compression 
and demonstrates that confocal microscopy is a powerful tool to decouple the ge- 
ometric and stress-bearing elements of the anomalous scaling of the bulk modulus 
in compressed emulsions. 

We would like to acknowledge Matthieu Wyart for pointing out the amplitude 
of the corrections to scaling and Eric Vanden-Eijnden for enlightening discussions. 
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